Quantum and quasi-classical types of superfluid turbulence 
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By injecting negative ions in superfluid *He in the zero-temperature limit (T < 0.5 K), we 
generated tangles of quantized vortex line with negligible large-scale flow. For this quantum regime 
of superfluid turbulence, the vortex line length L was found to decay at late time t as L ex. t~^; the 
prefactor being independent of the initial value of L. The corresponding effective kinematic viscosity 
is 0.1 K, where k is the circulation quantum. At T > 0.7 K, a jet of ions generates quasi- classical 
tangles identical to those produced by mechanical means. 

PACS numbers: 67.25.dk, 67.25.D-, 47.27.Gs 



In turbulent superfluids, flow on short length scales is 
restricted to quantized vortex lines, each of circulation 
K = h/m around a narrow core of radius ag, forming a 
dynamic tangle of the total vortex line length per unit 
volume L [J| . The turbulence can take two very different 
forms depending on whether the forcing is at scales above 
or below the mean inter- vortex distance (. = L~^l'^ . For 
flow on scales > £, the large quasi-classical eddies are 
the result of correlations in polarization of vortex lines 
0. On the other hand, when forced on quantum scales 
< the resulting uncorrelated tangle has no classical 
analogs and should have completely different dynamics 
first described by Vinen Q . In both cases, the dissipation 
of flow energy is through the motion of vortex lines; its 
rate per unit mass being |4| 
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However, the efficiency of the process, expressed through 
the "effective kinematic viscosity" v, can be different for 
these two regimes f^, which, following 0], we call "Kol- 
mogorov" (z^k) and "Vinen" (z/y) turbulences. 

In superfluid ^He, at high temperatures T > 1 K, 
where scattering of thermal excitations ("mutual fric- 
tion" ) dissipates the energy of the tangle and damps the 
waves on individual vortex lines (Kelvin waves), the val- 
ues of t'K and v\ for both types of turbulence are about 
the same, ~ 0.1k [1, HiQ- However, at sufficiently low 
temperatures where the dissipation can only be at very 
short wavelengths <C f [1] j to which the energy can be de- 
livered from larger scales ~ £ by a cascade of non-linear 
Kelvin waves [7[, ^vC?^) was so far unknown but vy^ was 
found to drop to 0.003k Q. 

To measure the values of t'y and z^k, one can monitor 
the free decay of homogeneous tangles. In any tangle, 
the quantum energy associated with the quantized flow 
on length scales r < £ is Eq — jL/ps (per unit mass) 0, 
where the energy of vortex line per unit length is 7 = 
BpsK^ and B w ln(£/ao)/47r is approximately constant. 
If the total energy is mainly determined by Eq, from 
Eq. Il]) we arrive at the late-time decay 

L = Bi^-H-\ (2) 

The typical values of £ are 0.1-1 mm. For ^He, K4 = 
h/nii = 1.00 X 10"'^ cm^/s and oq ^ 0.1 nm [3] mak- 



ing B « 1.2, while for ^He-B, K3 = /i/2m3 = 6.6 x 
10""* cmVs, oo ^ 13-65 nm [10] and B « 0.7. 

Now suppose that the tangle is structured due to the 
presence of flow on classical scales r > £, and the addi- 
tional energy of this classical flow Ec is much greater than 
-Eq. For the Kolmogorov spectrum between wavenum- 
bers ki and ^2 (^i *C ^2), while the size of the energy- 
containing eddy stays equal to the size of container d (i. e. 
ki « 27r/c?), the late-time decay becomes (6l. [Tl| 



L = {;iCf'^K-^k^^u^^'h-^/\ 



(3) 



The experimental value for the Kolmogorov constant for 
classical homogeneous isotropic turbulence is C « 0.5 x 
fl = 1.5(±11%) [12]. AH absolute values of j^k quoted 
in this paper are obtained from the fit to Eq. ^ using 
C = 1.5 and fci = 27r/d. 

In intermediate cases, the initial transient (before L(t) 
takes the form of either Eq. ([2]) or Eq. depending on 
initial conditions) can look much like L ccf^ with e « — 1 
0. However, the absolute value of the prefactor in the 
fit L cx t^^ will depend on the initial level of pumping 
and will never drop below Bv^^ . 




FIG. 1: (color online) Cartoon of the vortex configurations 
(side view) at different stages. The shaded areas indicate the 
trajectories of ions used to probe the tangle along two or- 
thogonal directions. To avoid contamination of the develop- 
ing tangle with new ions and vortices, only one probe pulse of 
ions was flred for each realization of the tangle, (a) At t = s, 
a pulse of CVRs is injected from the left injector. While most 
make it to the collector as a sharp pulse, some got entangled 
near the injector, (b) At t ~ 5 s, the tangle spreads into the 
middle of the cell, (c) At f ~ 20 s, the tangle has occupied 
all volume, (d) For up to 1000 s, the homogeneous tangle is 
decaying further. 

We generated turbulence in a cube-shaped container 
with sides d = 4.5 cm filled with ^He at p = 0.1 bar 
by injecting negative ions (electrons in a bubble state) 
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FIG. 2: (color online) Transients of the current of CVRs to 
the collector after different waiting times between the pulses 
(shown in seconds); the duration of injection is 0.3 s, V = 
90 V, T = 0.15 K. The arrival time (leading edge indicated 
by an arrow), 0.7 s, is the same for all pulses. The inset shows 
the dependence of the arrival time on the driving voltage. 



through the center of a square plate (Fig. 1). There were 
two field-emission injector tips ("left" and "bottom") and 
two collectors on opposite sides ( "right" and "top" ) [31 ■ 
The driving field of mean value V/d, shown in Fig. 1(a), 
was maintained by applying a potential — ^ to the in- 
jector plate and —V/2 to all side plates relative to the 
collector grid. The turbulence was detected by scatter- 
ing short time-resolved pulses of free ions at T > 0.7 K 
or charged vortex rings (CVR) at T < 0.7 K off the vor- 
tex tangle Q. The ion- vortex trapping diameter a for 
horizontal injection-collection was calibrated in situ on 
vertical vortex arrays of known density at steady rota- 
tion of the cryostat [8]. 

At T < 0.7 K, each injected electron dresses itself in 
a quantized vortex ring, somewhere between the tip and 
the nearby grid, and then propagates along with it into 
the cell. The ring energy Er — ^K^psR{Ai- — 2) and ve- 
locity Vr ~ j!^{A,- — 1) depend solely on its radius R 
(here — ln(8i?/ao) « 11) 0. The arrival of pulses 
of current carried by ballistic CVRs is shown in Fig. 2. 
The extrapolation of the dependence of the time-of-fiight 
on the driving voltage to V = gives 0.3 s (inset) corre- 
sponding to CVRs injected through the grid with initial 
velocities vq « 15 cm/s, radius Rq ~ 0.53 /zm and en- 
ergy 21 eV. Even though the side and bottom tips had 
very different threshold voltages for electron injection, 
300 and 120 V, these initial radii of the rings were about 
the same for both. No dependence of the ring energy on 
the injected current in the range 10~^^ - 10"^" A was 
found. The mutual friction a limits the range of a ballis- 
tic ring to i?o/Q: @ . For CVRs to survive the distance to 
collector d one needs a < R^/d ^ 10~^ (corresponding 
to T < 0.5 K), although a propeUing force due to the 
driving field extends this temperature range. 

When such rings collide, they build a tangle which then 
spreads into all volume. The decrease in the peak height 
/ of the current pulses due to the ballistic CVRs (shown 
in Fig. 2), relative to the long time amplitude I{oo) when 
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FIG. 3: (color online) Free decay of a tangle at T = 0.15 K. 
The injection direction and duration, and driving field are 
indicated. Probing with pulses of CVRs of duration 0.1-0.3 s 
were done in the same field as the initial injection, and also in 
the same direction except in one case (*). The line L oc 
corresponds to Eq. Q with za/ = 0.1k. 



all vortices have decayed, is the measure of blocking the 
paths of some ballistic CVRs with the probability per 
unit length La where a ^ R. We hence determine the 
vortex density after time t since stopping injection as 
L{t) — (crd)"^ In -j^. The slowly decaying component 
of the current (i.e. of the decay time much longer than 
the pulse width 0.1-0.3 s and the time constant of the 
amplifier 0.15 s), that appears in Fig. 2 as "tails" after 
each pulse at times t > 1.5 s, is due to the ions trapped 
on entangled vortices originating after each pulse. Their 
contribution, clearly seen in Fig. 2 as being of ~ 0.04 pA 
after only 5 s of waiting between pulses, was subtracted 
when determining I{t). The free decay for these tangles 
L{t) at T < 0.5 K is shown in Fig. 3, which is our central 
result. We found that the late-time decay is insensitive to 
the details of generating the initial tangle such as injec- 
tion duration (0.1 - Is) and intensity (10^^^ - lO'^" A), 
and driving field (0 - 20 V/cm) (see Fig. 3), i.e. the 
decay curve is universal and follows Eq. This gives 
strong support to the interpretation that the dynamics 
is that of the random tangle and not just the transient 
of a structured tangle with a not yet saturated energy- 
containing length. No changes with temperature were ob- 
served at T < 0.5 K. The values of ^v/k, obtained using 
Eq. ([2]), were slightly different for L{t) probed in horizon- 
tal (0.083 ±0.004) and vertical (0.120 ±0.013) directions. 
This might reflect the fact that trapping diameters a were 
slightly different for CVRs injected by different tips. As 
we could only calibrate a for the horizontal direction, 
and then used the obtained value to quantify L for both 
directions, the absolute value of vy — 0.08k seems more 
reliable. Note that, because of the finite width of the 
probe pulses 0.1-1.0 s and recorded time of flight ^^0.15 
- 1.0 s for CVRs and free ions, the decay time t has 
an uncertainty of 0.1-1 s; this is negligible for points at 
t > 20 s that were used for quantitative analysis. 
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During steady injection, the slow component of the col- 
lector current, associated with the spread of the tangle 
containing the trapped ions, only arrived after ~ 10 s. 
This implies that the initial tangle is always created near 
the injector and not near the collector at the opposite 
side of the cell. To study the dynamics of how the tan- 
gle spreads out into space, we probed L{t) across the 
direction of injection (Fig. 3(*)): within some 5-10 sec- 
onds the tangle reaches the center of the cell, then fills all 
volume and becomes homogeneous (as L{t) becomes in- 
distinguishable from those measured along the direction 
of injection) after ~ 20 s. The dynamics of spreading 
was found to be independent of the driving field. The 
observed time, ~ 20 s, is surprisingly short if one com- 
pares it with the estimate cP/D ^ 2 x 10^ s based on the 
simulated value of the diffusion coefhcient .iJl D ~ 0.1/v 
for L. However, if we consider "evaporation of vortex 
loops" from the boundary regions of strongly inhomoge- 
neous tangles, the rate of spreadin g is comparable with 
^ 1 mm/s observed in simulations [15|. Also, the initial 
tangle might maintain a certain polarization as a mem- 
ory of the orientation of the initial injected CVRs that 
it was created from. This would speed up the process of 
spreading before the tangle occupies all volume. 

The field due to the trapped charge was found not 
to affect the dynamics of the tangle's decay. This was 
checked for various ratios of the total injected charge of 
density n to L up to n/ L ^ 10^ cm~^, as well as in dif- 
ferent driving fields between —5 V/cm and 20 V/cm. As 
the time of flight of CVRs (Fig. 2) was independent of 
the trapped charge (but dependent on the driving volt- 
age) the field due to the trapped charge was always much 
smaller than the driving one. 




FIG. 4: (color online) Free decay of a tangle produced by 
a jet of ions (•) from the bottom injector into a 10 V/cm 
field for 150 s, as well as by an impulsive spin-down to rest 
d from 1.5 rad/s and 0.5 rad/s, at T = 1.60 K. All tangles 
were probed by pulses of free ions in the horizontal direction 
(spin-down data: probe field 20 V/cm, pulse length 0.5 s; ion 
jet data: probe field 10 V/cm, pulse length 1.0 s ). The line 
L oc t~^^^ corresponds to Eq. (|3]) with i^k = 0.2k. 

Yet, after a sufficient injection time the classical energy 



Ec can become dominant. Indeed, the observed late-time 
decay at T > 0.7 K (Fig. 4), L oc especially promi- 

nent for dense initial tangles after long (> 10 s) injection, 
was virtually identical to that of turbulence produced by 
mechanical means. Using Eq. ([3]), we extract the values 
of t'K(2^) for L oc t^^^^ and plot them in Fig. 5 for ion-jet, 
spin-down and grid turbulences The spin-down data 
include those published [8] and new points. The revised 
fitting routine, in which we only fit the L{t) for t greater 
than at least 50 s (to ensure that the transient behavior 
is not mistaken for the late-time decay), has resulted in 
slightly smaller the values of vk at T near 0.8-1.0 K than 
those published in Q . 
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FIG. 5: (color online) The effective kinematic viscosity for 
different types of superfiuid turbulence. Vinen (quantum) 
type, generated by: CVRs (this work, L{t) measured along 
the vertical (A) and horizontal (v) directions), build-up of 
counterflow 0], simulations [l^, vibrating grid in ^He-B vs. 
temperatures at which ''He has comparable values of 10~* < 
a < 2 X 10"^ [loj. Kolmo gorov (quasi-classical) type: ion 
jet (this work), spin-down sampled by CVRs (A) and free 
ions (v); towed grid pd| ]. 

Let us now discuss the temperature dependence of 
iA/{T) (Fig. 5). For T = 0, a tangle of vortex fila- 
ments initiated at short lengths and then allowed to de- 
cay through frequent reconnections of vortices was sim- 
ulated using the local induction approximation [l^ , and 
the L oc decay was observed yielding i/y = 0.06k. 
Another theoretical simulation, of the decay of a vortex 
ring through sound emission by Kelvin waves without 
any reconnections gives only i/y = 1.1 x 10~^k [iTl] thus 
emphasizing the importance of reconnections. Experi- 
mentally, only inhomogeneous tangles were generated in 
the T = limit by colliding vortex ring s before this work: 
by a jet of *He through an orifice [1^, injection of elec- 
trons ^9*1 and by vibrating objects such as a grid in ^He-B 
[20| . Only the latter was studied quantitatively and re- 
vealed that within a range of levels of initial pumping, the 
free-decay transients (before switching to late-time decay 
L oc t~^^'^) of type L (x were observed with a pref- 
actor increasing with increasing initial drive. Applying 
Eq. (21) to the lowest one hence yields the lower margin 
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oi lA/ = 0.27k3 for ^He-B which is factor 2-3 larger than 
our i/y for ''He. Of corse, the dissipation in ■^Hc-B might 
have additional channels [H- 

At T > 1 K (a > 10~^), computer simulations of tan- 
gle's decay while enforcing Vn = [H] yield an increase in 
vy{T) with temperature. The experimental vw, obtained 
for forced counterflow in narrow channels [s^ , are system- 
atically lower than the simulated ones. With increasing 
temperature the normal component becomes less viscous, 
hence it becomes impossible to enforce the = const 
condition. Yet, approaching T ^ 1 K from above (as the 
normal component firstly becomes laminar and eventu- 
ally irrelevant) the high-temperature data for seem to 
converge to the T — Q value ^ 0.1k. 

There exist two models for the drop in vk upon ap- 
proaching the r = limit. L'vov et al. ^] argue that 
conversion of the quasi-classical energy into Kelvin waves 
has constraints which cause a pile- up of vorticity at scales 
> i and hence a decrease in V]^. As Vinen tangles do not 
possess any energy at large scales, this scenario is consis- 
tent with vw{T) having no similar drop. Alternatively, 
Kozik and Svistunov [2l| explained fx {T) for T < 1 K as 
the result of fractalization of the tangle at smaller length 
scales < £ as a{T) decreases between ^ 10~^ and 10^^. 



To interpret the observed weak temperature dependence 
of vv{T) within this approach one has to assume that, 
for Vinen tangles, fractalization does not add a signifi- 
cant extra contribution to L for T < 1 K. 

To conclude, we generated a homogeneous quantum 
(Vinen) vortex tangle in the T = limit and found that, 
independently of the level of initial forcing, its late-time 
decay follows the law L oc with a universal pref- 
actor within more than two orders of magnitude. We 
hence extracted the effective kinematic viscosity for such 
a tangle, 0.08k, in reasonable agreement with the value 
0.06k from simulations of the reconnections-driven decay 
[l^, but about 30 times larger than that for a homo- 
geneous quasi-classical (Kolmogorov) turbulence (gener- 
ated by either impulsive spin-down to rest [sl or by ion 
jet in this work). Thus, in the zero temperature limit the 
energy flux fed at quantum length scales < £ generates 
little excess vortex length unlike the one fed at classical 
scales > £. 
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